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Structure of turbulent non-premixed 
flames modeled with two-step chemistry 

By J. H. Chen 1 , S. Mahalingam 2 , I. K. Puri 3 AND L. Vervisch 4 


Direct numerical simulations of turbulent diffusion flames modeled with finite- 
rate, two-step chemistry, A + B -* I, A + I -> P , were carried out. A detailed 
analysis of the turbulent flame structure reveals the complex nature of the penetra- 
tion of various reactive species across two reliction zones in mixture fraction space. 
Due to this two-zone structure, these flames were found to be robust, resisting ex- 
tinction over the parameter ranges investigated. As in single-step computations, 
mixture fraction dissipation rate and the mixture fraction were found to be statisti- 
cally correlated. Simulations involving unequal molecular diffusivities suggest that 
the small scale mixing process and, hence, the turbulent flame structure is sensitive 
to the Schmidt number. 


1. Introduction 

The development of turbulent combustion models that accurately reflect key phys- 
ical phenomena is essential for many engineering applications. Modern laser-based 
diagnostics applied to simplified laboratory flows are providing valuable data that 
is essential for any model development. However, since chemical reactions in flames 
proceed by a series of elementary steps, it is extremely difficult to isolate the role 
played by individual species and individual reaction steps on turbulent flame struc- 
ture. Thus the task of studying the two-way coupling between turbulence and 
combustion becomes formidable. Recently, direct numerical simulations (DNS) of 
turbulent non-premixed flames with simplified chemistry have proven to be useful 
in studying various aspects of the problem. 

DNS of turbulent diffusion flames including complex geometry and full chemistry 
is neither feasible with present-day computers, nor is it desirable since the objectives 
of DNS are usually to study specific aspects of the full problem. It is in this spirit 
that we studied the influence of turbulence on the structure of non-premixed flames 
in which chemistry is modeled through the following two- step mechanism: 

A + B^I ; A + /-P. (1) 

The stoichiometric coefficients are chosen so that the global step is identical to 
the one step mechanism discussed by Chen et al (1992). The mechanism and the 
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parameters were chosen to model stepl as a radical production step in which the 
intermediate I is produced followed by a radical consumption step in which the 
intermediate I is consumed to form the product P . Step2 proceeds with a small 
activation energy and a large enthalpy of reaction compared with stepl. Usually, 
the intermediate radical species I diffuses at a faster rate compared with the other 
species. Only Fickian diffusion is modeled. Such reaction mechanisms have been 
studied using large activation energy (for both steps) asymptotics by Margolis and 
Matkowsky (1982). They point out that the case wherein the radical species con- 
centration has a non-zero peak is typical of multi-step flames. The selected reaction 
mechanism may also be used to model practical flames (such as hydrocarbon oxida- 
tion) using global parameters. Altering the molecular diffusivity of the intermediate 
species can lead to modifications to flame structure and a change in flame temper- 
ature. This latter effect could be significant in accurately predicting thermal NOx 
formation in practical combustors (Law and Chung, 1984). Experimental mea- 
surements by Kerstein et al (1989) in nonreacting flows suggest that effects due 
to differential diffusion are not insignificant, and, hence, it is important to include 
them in models of chemically reacting flows. In reacting flows with significant heat 
release, these effects are likely to be amplified because of significant decrease in local 
Reynolds numbers due to an increase in kinematic viscosity with temperature. 

The goal of this work is thus to obtain a good understanding of turbulent non- 
premixed flames modeled by a two-step mechanism, to identify the role of the 
intermediate species on flame structure, and to understand the significance of dif- 
ferential diffusion on flame structure. 

2. Model problem 

The problem studied is the interaction of turbulence with an initially unstrained 
laminar diffusion flame. Computations to date have been performed over a two- 
dimensional square domain with 129 x 129 equi-spaced grid points. A laminar 
diffusion flame in which species A and B are present on either side of the flame 
(located at approximately x = 0) is initialized. Non-reflecting boundary conditions 
in x and periodic boundary conditions in y are imposed. The full compressible form 
of the equations of continuity, Navier-Stokes, energy, and species are solved using 
higher-order finite differencing schemes. The code used was developed by Ttouve 
(1991) for 3D premixed flames and subsequently modified for non-premixed flames. 
The treatment of boundary conditions and the differencing scheme is based on the 
methodology of Poinsot and Lele (1991). Fluid viscosity varies with temperature 
according to a power law with an exponent of 0.76. 

2.1 Governing equations and parameters 

The governing equations and non-dimensionalization are discussed by Chen et 
al (1992). In this subsection, the evolution equations of various species and the 
parameters associated with the two step chemistry model are discussed. 

The generalized chemical scheme is given by: 


v Ax A + vqB — ► uj x I ; ua 7 A + i// 2 Z — ► v P P, 


( 2 ) 
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where the u's axe the stoichiometric coefficients. 
We define the following mass ratios: 
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where W A , Wg, Wj, and W P are the molecular weights of species A, B, J, and P 
respectively. 

In terms of mass, the overall reaction is given by 


1A + rB — > (1 + r)B 


Wbvb 

W A [u Al +UA t j^\ 


(4) 


Let qi and q 2 be the energy released per unit mass of species A consumed in 
stepl and step2 respectively, and let 6 = 92/^1 • Then from the definition of the 
adiabatic flame temperature, one can show that the following equality, expressing 
global energy conservation, has to be satisfied: 


Qi [ e + £(1 ~~ e )) — TP ^ ^ c p(Tf ?o) 1 e “ , *2l ’ ^ 

*.4,0 V A\ ' 

where T f is the adiabatic flame temperature, Y Afi is the mass fraction of species 
A in the unmixed A stream, c p is the specific heat, To is the temperature of the 
unmixed A and B streams, and is the overall equivalence ratio defined as: 


<f> = r 


Y a , q 
V fl.o’ 


( 6 ) 


where Yg , o is the mass fraction of species B in the unmixed B stream. Continuity 
and Navier-Stokes equations are given in Chen et al (1992) and are not repeated 
here. The energy and species equations are: 
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where 


_ pUkUk 
2 



( 12 ) 


In the above equations, p denotes the fluid density, E denotes the total (internal 
plus kinetic) energy per unit mass, and 7 is the specific heat ratio taken to be 1 . 4 . 
The thermal conductivity A and the molecular diffusivities and V? 

are related to the dynamic viscosity p through the Prandtl and species Schmidt 
numbers. The other symbols have the usual meaning as in the standard literature. 
The reaction rate terms wa x and wa 2 are given by the Arrhenius expressions: 


UU, = BipY A pY B exp (~^r j , and w Aj = B 2 pY A pYiexp , (13) 

or alternatively in the form: 

w Al = Bj pY A pY B exp ( ~_^ a ( ( 1 1 ~ 7 ' T ) ) ) ’ (14) 

and 

w Al = B 2 pY A pY f exp ( t , (15) 

where T ai and T a2 are the activation temperatures associated with stepl and step 2 , 
r is a nondimensional temperature defined as r = (T — 7o)/(7/ — To) as given 
by Williams (1985), Bj = B\ exp(-/?i/a) and B 2 = B 2 exp(-/? 2 /a) are the pre- 
exponential factors, and a is a temperature factor defined as a = (T/ — Tq)/T/. 

In the simulations reported here, the following values were chosen for the stoichio- 
metric coefficients: va x = va 2 = 1 / 2 , vq = i //* = vj 2 = up = 1 , so that the overall 
scheme corresponds to the one-step model used by Chen et al (1992). Values for 
other constants are: a = 0.8, /?i = 8 , /? 2 = 2, 6 = 5, consistent with the discussion 
in Section 1 . The value of the Prandtl number is unity. The Schmidt numbers are 
Sca = Scb = Sep = 1 . The Schmidt number of the intermediate species Scj = 1 
unless otherwise noted. 

2.2 Initial conditions and turbulence field 

The generation of initial conditions corresponding to profiles associated with a 
one-dimensional, unstrained laminar diffusion flame was achieved by the following 
procedure: the flowfield is initialized with the solution to the one-step chemistry 
problem discussed by Chen et al (1992). The second reaction is turned “off” by 
setting B 2 = 0 , and the product profile for the one-step model is used as the 
profile for species I. The equations are time- advanced, and, simultaneously, the 
pre-exponential coefficient B 2 is increased until the desired value is reached. Care 
is taken to ensure that the flame does not extinguish during this process. Once the 
pressure waves exit the domain, the flowfield is saved (after any required rescaling), 
and, from this point forward, the turbulence is allowed to interact with the flame. 
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FIGURE 1. Species concentration profiles across the initial ID laminar flame for a 
representative condition (Da = 1). 

For details on the turbulence characteristics and definition of the global Damkohler 
number, Chen et al (1992) may be consulted. A plot of the species profiles of 
the laminar flame is presented in Figure 1. It is clear that the concentration of 
the intermediate is significant. Its peak value is a measure of the extent of phys- 
ical separation between reaction zones associated with the production (stepl) and 
consumption (step2) of the intermediate species. 

3. Results and discussion 

A quantity that is often used in descriptions of non-premixed flames is the mixture 
fraction, defined for our choice of stoichiometric coefficients, molecular weighty and 
unmixed stream concentrations as: 

z _ 2Y A -Y,-2Y B + 2 (16) 

4 

such that in the unmixed A stream, Z = 1, and in the unmixed B stream, Z = 0. 
When all the molecular diffusion coefficients are equal, Z is a conserved scalar. 
In the limit of infinitely fast chemistry, at the flame Y A = Y B = >7 = 0, giving 
Z ci = 1/2, whereas if the second step is relatively slow compared to the first step, 
at the flame Y A = Y 0 = 0, Y, = 1, giving Z c2 = 1/4. Thus these two limits of 
Z identify the locations of the reaction zones under these limiting cases. Finite 
rate chemistry will tend to cause a deviation from these two stoichiometric surface 
locations. Note that Z is not a conserved scalar when the diffusivities of the various 
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FIGURE 2(a). Distribution of species A mass fraction Y A , versus mixture fraction 
Z {Da = 1). 


species are not equal. In Figure 2(a), the mass fraction Y A is plotted as a function 
of Z. It is clear that most of the data points lie close to the equilibrium line, 
represented by the lower boundary of the data points. The penetration of species A 
into the B regime is due to finite rate kinetics. The scatter in these data points is 
due to unsteady effects associated with the turbulence field. The variations of Y B 
and Yj and the reduced temperature r are shown in Figures 2(b)-(d). These results 
are consistent with Figure 2(a). It is apparent from Figure 2(d) that, no significant 
drop in temperature is seen, suggesting that extinction is not seen for the conditions 
of the present simulation. 

Reaction rate profiles associated with stepsl and 2, respectively, are presented in 
Figures 3(a)-(b). Several features may be readily observed. It is clear that reaction 
rates for steps 1 and 2 peak for Z ~ 0.45 and Z ~ 0.58, respectively. Since the 
mass fraction of species I is significant where reaction 1 occurs, it is apparent from 
the definition of Z that the shift would be towards lower values of Z. Reaction 
2 occurs on the A side of the flame, where A is in excess, thus shifting the peak 
towards larger values of Z . For the parameters chosen, we have two reaction zones 
that are not completely segregated in mixture fraction (and physical) space. It is 
clear from these plots and the plot of reduced temperature (Figure 2(d)) that no 
extinction has occurred in either step. This is a significant observation since the 
global Damkohler number is comparable to the case for single step chemistry for 
which extinction was observed (see Chen ct al , 1992). An increase in the number 
of radical-like species and segregation of reaction zones apparently makes the flame 
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FIGURE 2(b). Distribution of species B mass fraction Yq , versus mixture fraction 
(Da = 1). 


0.25 



Mixture Fraction 

FIGURE 2(c). Distribution of intermediate species I mass fraction 1/, versus 
mixture fraction (Da = 1). 
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FIGURE 2(d). Distribution of reduced temperature r versus mixture fraction 
(Da = 1). 

less sensitive to extinction either by strain or flame shortening. 

A quantity that may be used to characterize reaction-diffusion zones is the scalar 
dissipation rate of the appropriate scalar, defined for species A as: 


Xa = 2V a VY a • VY a (17) 

In Figures 4(a)-(c), the scalar dissipation rates of scalars A , B, and the interme- 
diate species I are plotted as functions of mixture fraction. Results for species A 
are very similar to that for species A in the single step case (compare with Figure 
2(b) in Chen et al, 1992), with the peak occurring at a mixture fraction of approx- 
imately 0.7. Note, however, that there are no data points corresponding to frozen 
flow in the present case. This is consistent with the fact that no extinction was 
observed for the conditions simulated. Penetration of species A into the B side 
is apparent, but this is solely due to finite-rate kinetics. The scatter in the data 
points is due to unsteady effects. Results for species B complement those for A. 
Results for species I are interesting. If the second reaction were suppressed, it is 
clear that one would see a minimum in xi for Z = 0.5, with two symmetric peaks 
on either side, corresponding to locations where the molecular diffusion of species 
I would be highest. Note that since I is consumed by A in step2, the peak on the 
A side is enhanced, relative to the peak on the B side, with the minimum shifted 
to approximately Z = 0.4. This location corresponds to generation of species I by 
step 1. 
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FIGURE 4(c). Distribution of dissipation rate of intermediate species I versus 

s mixture fraction (Da = 1). 

i 

i 

In Figure 5(a), the distribution of the scalar dissipation rate of the intermediate 
as a function of the mixture fraction for Scj = 0.5 is shown. It is clear that the 
rapidly diffusing intermediate species penetrates further into the A side (compare 
i with Figure 4(c)). Instantaneous contour plots of in physical space for Scf = 1 

- and Scj = 0.5, shown in Figures 5(b)-(c), reveals that the dissipation zone is larger 

| when the species diffusion coefficient is larger (ie., Sc/ lower). This result suggests 

i that small scale mixing is sensitive to the molecular diffusion coefficient of the 

intermediate species even though turbulent diffusion is the more active process. It 
, is clear that three dimensional results need to be examined before one can draw 

5 definitive conclusions. 

' In Figure 6, the mixture fraction dissipation rate x defined as: 

X = 2WZ • VZ, (IS) 

is plotted as a function of mixture fraction. 

This picture is qualitatively similar to the results for single step chemistry pre- 
sented by Chen et al (1992), with the exception that no extinction is observed in the 

- present case. It is clear that a correlation exists between \ and % contrary to the 

= assumptions made in laminar flamelet theories (Peters, 1986, Warnatz and Rogg, 

1986). 

_ 4. Conclusions 

From the first study of non-premixed flames modeled by a two-step reaction 
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mechanism in a compressible, turbulent flow, two-dimensional databases were gen- 
erated. Finite values of the intermediate species mass fraction arise as a result of 
a slight segregation of the reaction zones associated with the two steps. A con- 
sequence of the segregated reaction zones is that the flame is less susceptible to 
extinction compared to flames modeled with a single-step, having an equivalent 
overall reaction step. It was demonstrated that the assumption of statistical in- 
dependence between the mixture fraction and its dissipation rate is a poor one, 
invalidating an important assumption made in flamelet models. It is expected that 
these results would hold for three-dimensioned simulations currently in progress. 
Preliminary results from simulations in which the intermediate species diffuses at 
twice the rate compared to the other species suggest that the small scale mixing 
process and, hence, the turbulent flame structure is significantly influenced. Fur- 
ther examinations of the flame structure and flame response to turbulence are being 
carried out. 
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